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In multifragmentation of hot nuclear matter, properties of fragments embedded in a soup of 
nucleonic gas and other fragments should be modified as compared with isolated nuclei. Such 
modifications are studied within a simple model where only nucleons and one kind of heavy nuclei 
are considered. The interaction between different species is described with a momentum-dependent 
two-body potential whose parameters are fitted to reproduce properties of cold isolated nuclei. The 
internal energy of heavy fragments is parametrized according to a liquid-drop model with density and 
temperature dependent parameters. Calculations are carried out for several subnuclear densities and 
moderate temperatures, for isospin-symmetric and asymmetric systems. We find that the fragments 
get stretched due to interactions with the medium and their binding energies decrease with increasing 
temperature and density of nuclear matter. 
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I. INTRODUCTION 

In energetic nuclear collisions, the participating hot 
nuclear matter after an initial dynamic stage of compres- 
sion expands to a subsaturation density and then disas- 
sembles into many fragments due to growing instability. 
Statistical models of different genres (T-Q have generally 
been successful in explaining the many features associ- 
ated with the fragment multiplicities, the caloric curve, 
the density of the fragmenting systems, etc. They also 
offer a broad hint about the general nature of the phase 
diagram of nuclear matter 0,111 at temperature T ~ 3— 8 
MeV at subsaturation densities p ~ 1/20 to l/5th of 
the normal nuclear density po- Possible liquid-gas phase 
transition and associated condensation to form nu- 

clear clusters at these temperatures and densities help in 
a better exploration of many phenomena of astrophysical 
interest like supernova explosions or explosive nucleosyn- 
thesis Hunt. 

Analysis of recent laboratory experiments [HI, [l9[ on 
nuclear multifragmentation seems to indicate that the 
properties of the nuclides arc modified at subnuclear den- 
sities (p ~ po/3) they are created in corresponding to 
freeze-out. The symmetry energy, for example, is re- 
ported to be progressively reduced [2(| with excitation 
energy, which is attributed to the in-medium modifica- 
tions of the properties of the hot fragments [U [22| . A 
looming uncertainty about whether the measured sym- 
metry energy corresponds to the hot fragments or the 



equivocal decision about the medium modifications. The 
reduction in symmetry energy can have a fair explana- 
tion from the thermal and expansion effects of the disas- 
sembling system [24 |. The surface properties of the hot 
fragments [2l| as well as their bulk energy [26| are also 
speculated to be modified due to the embedding envi- 
ronment. A quantum statistical approach to the nuclear 
equation of state taking into account the formation of 
clusters [U [H| shows that the properties of these clus- 
ters are modified due to the medium in which they are 
formed. The symmetry energy of low density warm nu- 
clear matter predicted by this model seems to be in good 
agreement with the experimental data pfj| . In a recent 
experiment _3Q] , it is claimed that the binding energies of 
very light clusters (A < 4) produced in multifragmenta- 
tion progressively tend to zero in the temperature range 
of T ~ 5 — 10 MeV even at a very low in- medium density 
- 0.05 po- 

The changes, if any, of the bulk properties of the frag- 
ments produced in nuclear disassembly are expected to 
originate from the effects of the residual interaction of 
the fragments with the surroundings. The aim of the 
present article is to study these effects starting from an 
effective nucleon-nuclcon interaction. To keep the physics 
simple and transparent and yet retaining all the basic es- 
sentials, we allow the dilute matter to condense into only 
one kind of nuclear species surrounded by a hot nucleonic 
gas and species of the same kind. Then we introduce the 
interactions between them and look for the minimum of 
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of the fragments maintaining chemical equilibrium be- 
tween the fragments and the nuclcon gas. The energy 
and free energy of the nuclei are evaluated in the liquid- 
drop framework that makes it easier to account for the 
associated changes in the surface and symmetry energy 
with the change in the volume of the fragment species. 

The organization of the paper is as follows. In Sec. II, 
the outlines of the theory are given. Results and discus- 
sions are contained in Sec. III. The concluding remarks 
are presented in Sec. IV. 

II. THEORETICAL FORMULATION 

As is well known, hot low-density nuclear matter 
condenses into nuclear fragments of different sizes sur- 
rounded by nucleons. We postulate that the nucleons 
and the fragments interact through a common effective 
interaction. The bulk properties of the fragments may be 
modified because of this interaction too. For a qualita- 
tive understanding of such system, we take only one kind 
of fragment species of mass A and charge Z. The interac- 
tion is chosen to be the modified Seyler-Blanchard (SBM) 
interaction. Its properties are summarized in Sec. IIA. 
In Sec. IIB, the method for evaluating the nucleon- 
fragment and the fragment-fragment interactions is de- 
scribed. In Sec. IIC, we study the observablcs sensitive 
to the medium modification of the properties of finite 
nuclei. 



TABLE I: The parameters of the effective interaction (in 
MeV fm units) 



Q 


C'u 


a 


b 


d 


K 


348.5 


829.7 


0.6251 


927.5 


0.879 


1/6 



and k are listed in Table I. The procedure for determining 
these parameters are given in detail in Ref. [33|. These 
parameters are somewhat different from those given in 
[32l ] ; in the latter, the symmetry coefficient a sym for infi- 
nite nuclear matter at saturation density po was chosen 
to be 34 MeV; in the present calculation, it is taken to 
be 31 MeV to be more consistent with the recent esti- 
mates [H, HH . The effective mass of the nucleon coming 
from the momentum dependence of this effective inter- 
action is 0.62m for symmetric nuclear matter, where m 
is the nucleon mass. For the interaction, the isoscalar 
volume incompressibility A^, symmetry incompressibil- 
ity K sym and L, a measure of the symmetry pressure are 
240, —101, and 59.8 MeV, respectively. It is interesting to 
note that the symmetry coefficients a sym , L, and K sym of 
this interaction are within the range of values suggested 
by the empirical constraints emerging from recent analy- 
sis of different observables [3^ - [4pj . With this interaction, 
for symmetric nuclear matter, the critical temperature is 
reached at T c =14.9 MeV, when the surface energy van- 
ishes. 



A. The effective interaction 

The SBM interaction is a momentum and density de- 
pendent effective interaction of finite range. In the con- 
text of the nuclear mass formula, an interaction of this 
type has been used with great success by Myers and 
Swiatecki [3l[. It also reproduces the rms radii, charge 
distributions, and giant monopole resonance energies for 
a host of even-even nuclei ranging from 16 O to very heavy 
systems [32l |. Its form is given by 



B. Nucleon- fragment and fragment- fragment 
interaction energy 

A low-density nucleonic matter of density pb and asym- 
metry X = (p b n — Pp)/pb breaks up to a system of free 
(unbound) neutrons and protons of density p n and p p and 
a collection of mass- A fragments (we call it AN matter), 
all at temperature T. The clusterized matter is thermo- 
dynamically more favorable than the uniform matter at 
the same (low) density, asymmetry and temperature [4lj | . 
The baryonic density is then given by 



v{r,p,p) = C itU [vi(r,p) +v 2 (r,p)] , 
P 2 

"l = -(1 - pj)/(ri,r 2 ), 

V2 = rf 2 [^i) + ^ 2 )] K /(ri,r 2 ), (1) 



with 



/(ri,r 2 ) = 



3 -|ri-r 2 |/a 



ri - r 2 \ a 



(2) 



The subscripts I and u to the interaction strength C refer 
to like-pair (nn or pp) and unlike-pair (np) interactions, 
respectively. The relative separation of the interacting 
nucleons is r = r x r 2 and the relative momentum is 



Pb = Pn + Ap A , 



(3) 



where pn = Pn + Pp is the free nucleonic density and pa 
is the number density of the fragment species. The total 
thermodynamic potential f2 of the system is 



O 



TS - ^ VrNr - fi A N A , 



(4) 



where E, S, p T , \ia, A T ,and Na are the total energy, en- 
tropy, chemical potentials of the free nucleons and the 
fragments, free nucleon number and the number of the 
fragments of mass A, respectively. The isospin index 
(n,p) is represented by r. Chemical equilibration ensures 
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where N and Z are the neutron and proton numbers in 
the fragment. The total internal energy of the AN system 
is written as 



E — Enn + Ean + Eaa- 



(6) 



In Eq. (6), Ean is the contribution coming from the 
nucleon- fragment interaction (Van)- Enn measures the 
kinetic energy of the free nucleons plus the interaction 
energy amongst themselves. Eaa is the sum total of the 
the kinetic energy of the fragments, interaction energy 
among them and their binding energies. 

Assuming for simplicity that the fragments are sharp- 
surface liquid drops with a uniform nucleon density pi, 
these terms can be explicitly written as, 



E 



NN 




P 2 

cZridpi— — n r (pi) 
2m T 

x[ui(|ri - r 2 |, |pi - pa|) - 



+ 7j J dr 1 dp 1 dr 2 dp 2 

v 2 (|ri - r 2 \,2p N )} 



x[C;n r (p 2 ) + C u n_ r (p 2 )]n r (p 1 ) 



(7) 
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FIG. 1: Space coordinates shown for nucleon (located at B) 
and fragment (with center at A) configuration. The origin of 
the coordinate system is at O and P is any arbitrary point 
within the fragment. 



Ean 



and, 



where, 



E AA 



\ J dr 1 dp 1 dr 2 dp 2 

xn T (Y 1 ,p 1 )h A {Y 2 ,p 2 ) S ^{Ci5 TTl +C U (1 -S TT i)) 

t' 

X lv dp ^'^' (r ' p ^' ) [ Wl(|r + R| ' |pi 



<pf T > + Pa)|) + v 2 (\r + R\,p N + pi) 



(8) 



(9) 



E AA = E AA -N A B A (pi,T), 

P 2 1 f 

dvdp- h A (r,p) + - / dridpidr 2 dp2 

2m A 2 J 

xn A {r ll pi)n A (r 2 ,p 2 ) / drdr' 
Jv A 



xnf T (r,pf T )hf T ,(r',p£,) Wl (|R + r -r'|, |( Pl 

+P^)-(P2+P^')I) 



-u 2 (|R + r-r'|,2^ 



(10) 



In Eq. (8), R = r 2 ri is the distance between the nu- 
cleon and the center of the nucleus (see Fig. 1). In 
Eq. (10), R is the distance between the two fragment 
centers (see Fig. 2). The various space coordinates oc- 




FIG. 2: Space coordinates shown for fragment-fragment con- 
figuration with O as the origin of the coordinate system. P 
and P' axe arbitrary points within the fragments with A and 
A' as their centers. 



respectively. In evaluating the coordinate-space integrals 
in Eqs. (7), (8) and (10), we have assumed, as in the cal- 
culation of the equation of state of dilute nuclear matter 
[42l HU that the free nucleons do not penetrate the sharp 
surface nuclei and also that the fragments do not inter- 
penetrate so that the identity of the free nucleons and 
the fragments is never altered. This 'no-overlap' approx- 
imation is tantamount to use of the 'excluded-volume' 
correction employed earlier [H where the 'free' volume 
available to fragments is reduced compared to the total 
volume V by at least the internal volume of the nucleons 
and the fragments. 

In the above equations, m T and m A are the masses 
of the nucleons and the fragments, and B A (pi,T) is the 
binding energy of the produced fragments at temperature 



4 



ft-A = |I«a; n r and riA are the occupation probabilities 
for the free nucleons and the fragments and qa is the de- 
generacy of the fragments taken to be 1 or 2 depending 
on whether they are bosons or fermions. Since the system 
is infinite, the occupation functions h T are independent 
of space coordinates. Similar is the case for tia- In the 
equations above, space dependence has, however, been 
retained to correlate with Figs. 1 and 2. Since the frag- 
ments are also taken to be uniform drops, nf T {— Tsn^.), 
the distribution function of the constituent nucleons in- 
side the fragments is also independent of the space coor- 
dinates. The explicit functional dependence of the distri- 
bution functions on the space coordinates is henceforth 
omitted from the equations where the distribution func- 
tions may enter. The momenta of these nucleons inside 
the fragment is designated by pf r . The notation j v 
refers to the configuration integral over the volume Va of 
the fragment. All other integrals are over the entire con- 
figuration or momentum space unless otherwise specified. 
The distribution functions yield the densities as, 



2 



J n T {p)dp = N T /V = pr, (11) 



J n A {p)dp = N A /V = p A , 



2 



J nf T (p)dp = A T /V A = plr, 



(12) 



(13) 



where V is the volume of the AN matter, A T is the neu- 
tron number N or proton number Z in the fragment, 
Pn = Y. t Pt, Pi = Y^rPir and V A = ^-R\- In Eq. (13), 
pi T refers to the neutron or proton number density in the 
fragment; Ra is its sharp-surface radius. 
The total entropy of the AN system is 



>J A 



oint 



(14) 



where in the Landau quasiparticle approximation, the 
contribution ^ T &t of the free nucleons is taken as 



h 3 ^ 



drdp 



n r (p)lnn T (p) 



+(l-n T (p))ln(l-n T (p)) 



(15) 



ST is the entropy from the center-of-mass motion of 



the thermodynamic potential can be absorbed alongwith 
the binding energy term of Eq. (9) in the free energy of 
the fragments F A (pi,T) = (-B A (pi,T) - TSf*) when 
the thermodynamic potential takes the form, 



ft — Enn + Ean + E° A a — T {^2, ^ T ^) 

r 

- PrN T - p A N A + N A F A (p u T) (17) 



Minimization of fi with respect to n T and tia, remem- 
bering that Sn T (p) and 5riA(p) are separately arbitrary 
over the whole phase space, after some algebraic manip- 
ulations, yields for the distribution functions with the 
following structures, 



n T {p) 



n A (p) 



exp 



exp 



2mZT 



+ 1 



2m* A T 



-Va))±1 



(18) 



(19) 



In Eqs. (18) and (19), rj T = [p T -V?-V?)/T and T] A = 
(pA — Fa — Va)/T are the fugacities pertaining to the free 
nucleons and fragments, respectively and m* and m* A are 
the effective masses of the nucleons and the fragments in 
the medium, the masses getting renormalized owing to 
the momentum dependence of the force. The nucleonic 
rearrangement potential V 2 originates from the density 
dependence of the interaction. The effective nucleon and 
fragment masses are given by 



and 



— + 2V T ' 



m A 



(20) 



(21) 



where p 2 V^ and p 2 V\ are the momentum dependent con- 
tributions to the single-particle potentials V T and Va ■ 



Vr( P ) = V°+p 2 V T 1 , 



Va(j>) 



(22) 



(23) 



the fragments, and ST is their internal entropy. SI is ^ * , . , , , , 

° A A Expressions tor tnc momentum independent components 



evaluated as 



-|j J drdp[n A {p) lnn^(p) 
±(lTnx(p))ln(lTnA(p))]. (16) 



In the above equation, the upper and lower signs cor- 
respond to a fermionic and bosonic fragment, respec- 



V® and , alongwith those for V^ , V? and V\ are given 
in the Appendix. 

C. Energy and free energy of the system 

We take recourse to liquid-drop model for the cvalu- 
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mass A, charge Z and neutron number N at a con- 
stant density p and temperature T. The energy Ea(p, T) 
(=-B A (p,T)) is given by 



E A ( P ,T) 



a v {p,T)A + a s { Pl T)47rR A A 2 / 3 
(N — Z) 2 

+a sym {p,T) 
-Z 2 e 



.4 



1 



V Ra Rws 



(24) 



The term a v is the volume energy term for symmetric 
nuclear matter. Alongwith a v , the surface energy coef- 
ficient a s and the symmetry energy coefficient a sym are 
all density and temperature dependent. The last term in 
Eq. (24) is the the Coulomb term. One may note that 
the Coulomb energy is different from that for an isolated 
nucleus. As the fragment is embedded in clusterized 
matter, its Coulomb energy gets 'dressed'. It is calcu- 
lated in the Wigner-Seitz approximation 0. Here Rws 
is the radius of the spherical Wigner-Seitz cell, given as 
Rws = (f ^Pa)" 1 / 3 - The Coulomb energy has no ex- 
plicit temperature dependence. The radius Ra of the 
liquid drop is given by R A = ^l 1/3 /[|7rp(T)] 1 / 3 . In a 
similar vein to Eq. (24) , the free energy of the nucleus is 
taken as 



F A (P,T) 



f v (p,T)A + f s (p,T)4irR A A 2 / 3 
(N - Zf 



+f sym (p,T) 
-Z 2 e 



A 



i 



V Ra Rws 



(25) 



The volume terms a v and f v are calculated for symmet- 
ric nuclear matter at density p and at temperature T 
employing the SBM interaction. The density and tem- 
perature dependence of the surface free energy coefficient 
is assumed to be factorized \4M and is taken as 



fs{p,T) =a s (p o ,T = 0)U(p)y(T), 



(26) 



where a s (po , T = 0) is the surface energy coefficient at 
nuclear matter saturation density po at T =0. The ex- 
pressions for U(p) and y(T) are taken from Refs. [45| 
and [H , respectively. They are given as 



U{p) = l-^( P 



Po 



Po 



and 



y(T) = 



(T 2 — T 2 \ 5 / 4 
\T 2 +T 2 



(27) 



(28) 



T c is the critical temperature for nuclear matter cal- 
culated to be 14.9 MeV with the SBM interaction. 
The value of a s (p a ,T = 0) and k p are taken to be 
1.15 MeVfm~ 2 and 5.0, respectively. The surface entropy 
per unit area S sur f is obtained from f s as 



5- 



df s 



(29) 



which yields 



a s (p, T) = f s (p, T) + TS sur f 

, T 2 T 2 



a s (p ,0) y(T) 
T 2 T 2 



T 2 + T 2 



1/4 



( T 2 +T 2y 



U{p). 



(30) 



The symmetry coefficient a^ ym is dependent on the 
nuclear mass. It is taken as 1351 



a sym {p = po,T = 0) = 



a 4-1/3 ' 



(31) 



where a is the symmetry coefficient of cold symmetric 
nuclear matter taken as 31.0 MeV and ^ =2.4. For infi- 
nite matter, it is generally seen that a syrn decreases with 
temperature whereas f sym shows the opposite tempera- 
ture dependence j4(|. A nearly similar trend has been 
observed for finite nuclei |47j |: here f sym increases with 
temperature (though in some cases, an occasional de- 
crease is seen at low T). The density dependence of the 
symmetry energy of nuclear matter calculated with the 
SBM interaction is seen to be given by ~ (p/po) 7 with 
7 ~ 0.69 24|, in consonance with the reported exper- 
imental behavior (20| . With this in mind, we write the 
symmetry free energy coefficient in a factorized form as 

fs V m{p,T) = a sym (p,T= 0)g(T), (32) 

where 

a sym {p,T = 0) = a sym {p ,T = 0)(p/p o ) J . (33) 
For g(T), we assume a polynomial in T of the form, 



g(T) = (l + Vl T + V2 T 2 + Vi T A ). 



Then, 



Of, 



s lj III 



sym 



dT 



and therefore 

ym(Pi^) — fsymxPi'-^) 



T 



dfsym(p,T) 



dT 



(Is 



t (p,T = 0)[1 - is 2 T 2 -3i/ 4 T 4 ]. 



(34) 



(35) 



(36) 



(37) 



In a schematic model [47| , the observed T- dependence of 
a sy m and f sym has been seen to be moderately explained 
with values of 1^1,^2, and v± as —0.00848, 0.00201 and 
0.0000147, respectively, the dimensions of these quatitics 
being in relevant inverse powers of MeV. 

The internal entropy S 1 ™* of the fragments has con- 
tributions from the volume, surface, and the asymmetry. 
The latter two contributions have already been taken into 

j_ j_l l_ /on\ 1 /or\ l__ m 
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the fragments are taken to have uniform density p, the 
volume entropy is calculated using the expression given 
in Eq. (15) for symmetric nuclear matter at temperature 
T and at a density p. 

Combining the terms given by Eqs. (7)-(10), the total 
energy of the (n,p,A) system can then be written as, 



E = V 



x (1 - m*V T *) + -p T V T ° + PA TC A 



x(l - m* A V\) + - PA Vl - PA B A (p u T) 



.(38) 



We may remind here that p T corresponds to the free nu- 
cleonic density after condensation. The total entropy of 
the free nucleonic matter is, from Eq. (15), 

^2s t = V^2p t [-J 3/2 (vt)/Ji/2(Vt)-Vt ■ (39) 



Similarly, the translational entropy from the fragments 
is, 



S t X = N A [-C A -r ]A ], 
where C A is given by, 

C.4 = J3/2(VA)/Jl/2{VA), 

or 



(40) 



(41) 
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FIG. 3: (Color online) The percentage of nucleons contained 
in the fragments at a given baryon density pt at X = 0.0 
shown as a function of temperature after fragmentation. The 
left panels correspond to pt = 0.005 fm~ 3 , the right panels 
to pt = 0.02 fm -3 . The fragment specimens chosen are 40 Ca, 
56 Fe and 150 Sm, respectively. The blue full line corresponds 
to the case (1,1,1), i.e., the calculation where all three in- 
teractions NN, AN, and AA are included. The red dashed 
line refers to the case (1,0,1) where the AN contribution is 
neglected, whereas the black dotted line corresponds to the 
(1,0,0) calculation without both the AN and AA contribu- 
tions. For details, see text. 



C A = B 3/2 (i lA )/B 1/2 (r] A ), (42) 

depending on whether the fragments are fermionic or 
bosonic. In the above equations, the quantities Jk and 
Bk are the Fermi and Bose integrals; their definitions are 
given in Eq. (A. 4) in the Appendix A. As the fragment 
densities are usually very low, C A ~3/2. Since the in- 
ternal entropy S*™* of the fragments is now known as ex- 
plained earlier, the free energy T of the total AN matter 
can be calculated. 



III. RESULTS AND DISCUSSIONS 

In this paper, our primary aim is to investigate the 
changes in the properties of nuclei embedded in a hot 
medium of nucleons and other fragments produced in nu- 
clear multifragmentation. To simplify the problem yet 
retaining the main physics essence, we assume that after 
nuclear disassembly, the system contains a collection of 
only one kind of fragments of mass A and charge Z in 
thermodynamic equilibrium, with a hot soup of neutrons 
and protons. To begin with, we take a baryon matter 



asymmetry X. The binding energies and the free en- 
ergies of the nuclear fragments that enter into the cal- 
culation have been modeled in the context of the liquid 
drop mass formula. For the effective interaction, the mo- 
mentum and density dependent SBM force as scripted in 
Eqs. (1) and (2) has been chosen. Assumptions are made 
that the free nucleons do not penetrate the sharp-surface 
nuclei and that the fragments do not overlap. 

The three unknowns in the calculation are the free nu- 
cleon densities p n , p p and the fragment densities p A in the 
matter. The three constraints are the conservation of the 
total baryon number, the total isospin and the condition 
of chemical equilibrium between the nucleon gas and the 
fragments. For a given set of Pb,T and X, the calcula- 
tions start with a chosen value of the density pi of the 
constituent nucleons in the nuclear fragments at temper- 
ature T and an input guess density of p A . The energies 
and the free energies of the fragments are then known 
from Eqs (24) and (25). Exploiting the constraints, the 
final densities p A ,p n and p p are determined iteratively. 
The total free energy T of the given AN matter (see 
Eqs (6) and (14)) is then calculated as outlined earlier. 
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FIG. 4: (Color online) The nucleon number density in the 
fragments shown as a function of temperature, for baryon 
density p b = 0.005 fm" 3 at X = 0.0. The blue full line, 
red dashed line and the black dotted line convey the same 
meaning as in Fig. 3. The dash-dot (black) lines refer to the 
nuclear density when the nuclei ( 40 Ca, 56 Fe or 150 Sm) are in 
phase equilibrium with their own vapor. 



determines pi at a given baryonic density pb with asym- 
metry X at temperature T, which then, in the liquid 
drop framework, determines all the properties under in- 
vestigation of the produced fragments. 

The calculations have been done in the temperature 
range 2.5—8 MeV for symmetric and asymmetric nuclear 
matter. Initially a symmetric matter (X = 0.0) of low 
baryon density pb = 0.005 fm -3 is chosen, and the calcu- 
lations are then repeated at a higher baryon density pi, 
=0.02 fm" 3 . This helps to see how a denser medium ac- 
centuates changes in the nuclear properties. Three repre- 
sentative fragments are selected, namely, 40 Ca, 56 Fe and 
a heavier one 150 Sm. In Fig. 3, the percentage of nucleons 
in the fragments produced (^4/9a/p&x100) is shown as a 
function of temperature for the three fragment species. 
The left panels correspond to pb =0.005 fm -3 , the right 
panels refer to the higher baryon density pb =0.02 fm -3 . 
The blue full line is obtained from calculations with inclu- 
sion of all three interaction contributions, namely, (NN), 
(AN), and (AA) [see Eqs. (7)-(10)]. We refer to these cal- 
culations as (1,1,1). The dashed red line corresponds to 
calculations without the (AN) contribution and the dot- 
ted black line is the one obtained when both (AN) and 
(AA) contributions are excluded. The latter two calcula- 
tions are referred to (1,0,1) and (1,0,0), respectively. At 
low temperatures, most of the nucleons are contained in 
the fragments and the free nucleons are rare. This is ex- 
pected, vapor tends to condense to drops at low temper- 



probability decreases. Fragment formation also depends 
on the total baryon density pb\ at the higher pb, fragment 
formation probability is higher. As the system heats up, 
this probability goes down. 

Examination of the figure reveals few further features. 
The nucleus-nucleus interaction (AA) does not have a 
very significant role (as seen from the almost overlap- 
ping of the red dashed (1,0,1) and black dotted (1,0,0) 
lines), the nucleon-nucleus interaction (AN) is important, 
the importance grows with increasing baryon density pb- 
Normally, it is seen that the (1,0,1) or (1,0,0) calcula- 
tions favor the production of fragments compared to a 
full (1,1,1) calculation. These results allow us to conclude 
that heavy nuclei embedded in the medium arc affected 
mostly by the nucleons (and perhaps light clusters) sur- 
rounding them. This means that the description of the 
multicomponent nuclear system can be simplified by sub- 
dividing it into non-interacting cells containing one heavy 
nucleus and a proportional amount of the medium. This 
Wigner-Seitz approximation is widely used for studying 
inhomogeneous phases of nuclear matter [481450T ] . 

In Fig. 4, the progressive changes in the internal nu- 
cleonic density in the three nuclei produced from disas- 
sembly of symmetric nuclear matter of density ph=0.005 
fm -3 arc displayed as a function of temperature in the 
three panels. The blue (full line), red (dashed) and the 
black (dotted) lines have the same meaning as stated 
earlier. The dash-dot black line refers to calculations for 
isolated hot nuclei ( 40 Ca, 56 Fe or 150 Sm) in phase equi- 
librium with their own vapor (n-p gas). In that case, the 
phase-equilibrium conditions [5 1| determine the density 
and asymmetry of the embedding nuclconic gas along- 
with the internal density of the dipped nucleus. The 
dot-dash black line is drawn so as to serve as a reference 
against which the other ones can be compared. In this 
case thermodynamic equilibrium ensures that the tem- 
peratures of the nucleus and the surrounding (n,p) gas 
are the same, the pressure exerted by the nucleus bal- 
ances that of the gas and that the chemical potentials 
of the neutrons and protons of the nucleus are the same 
as those of the free neutrons and protons in the gas, re- 
spectively. Discussion on this part of p hysics is left out 
here, it is given in detail in Ref. (UHJ. The densities of 
the embedding n-p gas in clusterized matter and in the 
case of isolated hot nucleus are somewhat different. In 
both cases the density of the surrounding n-p gas is low 
at low temperatures and increases with increase in tem- 
perature. The asymmetry of this gas is also different in 
both cases. As an illustrative example, these properties 
of the embedding nucleonic gas are displayed in Fig. 5. 
The asymmetry of the disassembling system is X =0.2 
and the fragment concerned is 150 Sm. The asymmetry 
of the gas has a very insignificant role to play on the 
properties of the fragments as will be shown later. 

From Fig. 4, it transpires as expected, the fragment nu- 
clei in equilibrium in an embedding medium expand with 
temperature in all four calculations displayed in this fig- 
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FIG. 5: (Color online) The embedding nucleon gas density 
plotted as a function of temperature in AN matter where A 
is 150 Sm. The panels (a) and (b) correspond to pf, = 0.005 
fm -3 and 0.02 fm -3 , respectively. The red dotted, dashed 
and full lines refer to proton and neutron densities and their 
sum in the calculation for disassembled matter at asymmetry 
X —0.2. The black dotted, dashed and full lines correspond 
to proton and neutron densities and their sum in the phase- 
equilibrium calculation for 150 Sm. 



nuclei like 40 Ca or 56 Fe bloat up in volume by ~ 17- 
19% from their ground state equilibrium values when 
the temperature is raised to ~ 7—8 MeV. In the full 
(1,1,1) fragmentation calculation for nuclear matter at 
P6=0.005 fm~ 3 they do so by ~ 25 %. As is seen from 
the figure, incorporation of the nucleon-fragment (AN) 
interaction plays a significant role in the expansion of 
the fragments; also it is seen, as stated earlier, that the 
fragment-fragment interaction (AA) has little effect. In 
panel (c) of Fig. 4, it is noticed that the black dash-dot 
line does not extend beyond 6.1 MeV. This is the limiting 
temperature as obtained for 150 Sm in the phase equilib- 
rium calculation. Solutions for calculations without in- 
corporation of both AN and AA (1,0,0) interactions and 
for calculations without AN interactions (1,0,1) could not 
be obtained for this nucleus beyond 7.2 and 6.9 MeV, 
repectivcly. The so-found absence of solutions possibly 
points to the nuclear instability beyond these tempera- 
tures in these calculations. When the baryon density pb 
is increased, the nuclear soup of fragments and nucleons 
becomes denser. In this denser environment, a further 
expansion of the fragments, by another ^4% could be no- 
ticed. This is displayed in Fig. 6 where pb = 0.02 fm -3 . 
The (1,0,1) and (1,0,0) calculations for 150 Sm could now 
be extended to 8.0 MeV. 

As the internal density of the nuclei change progres- 
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FIG. 6: (Color online) Same as in Fig. 4 when pb = 0.02 
fm- 3 . 



reference calculation as mentioned earlier, it is expected 
that there should be a corresponding change in the bind- 
ing energies and free energies of the nuclei embedded in 
medium. We display those quantities for 40 Ca in Fig. 7. 
The left panels (a) and (b) show the results for the lower 
baryon density p&=0.005 fm~ 3 , the right panels (c) and 
(d) do so for p 6 =0.02 fm" 3 , both at X =0.0. The upper 
panels display the binding energy, the lower panels the 
free energy. The (1,0,0) calculations being nearly indis- 
tinguishable from the (1,0,1) calculations are not shown 
here. The results for the binding energies from the (1,1,1) 
and (1,0,1) calculations at the lower baryon density are 
nearly indistinguishable from each other, at the higher 
baryon density (ph=0.02 fm -3 ) only little changes are ob- 
served at the high temperatures. Compared to the refer- 
ence calculation (dash-dot black line representing phase 
equilibrium in a n-p gas), there is a gain in the bind- 
ing energy. This comes mostly from the decrease of the 
Coulomb energy because of the presence of other frag- 
ments (see Eq. 24). At higher baryon density, the frag- 
ment density is comparatively higher which explains the 
larger gap in binding energy with reference to the phase 
equilibrium calculation. For a given baryon density the 
reduction in the density of the fragments with tempera- 
ture reduces the gap. Similar arguments follow for the 
lower free energy of the embedded fragments as compared 
to that from the reference calculation. Not much of a dif- 
ference is seen for the heavier nucleus 150 Sm as displayed 
in Fig. 8. 

In order to see the importance of asymmetry on the 
fragment observables, the calculations have been re- 
peated for asymmetric nuclear matter with X = 0.2. In 
Fig. 9, comparison of the percentage of nucleons con- 



9 



10 



> 

<d 



< 

pq 



> 



< 



6 
4 
-10 
-12 
-14 
-16 



1 1 1 1 1 1 1 1 1 1 

\ 40 Ca 

" (a) 

- p b =0.005 (fm" 3 )\ 
1,1,1,1,1, 


I i I i I i I i I i 

" (c) s \\ - 

-p b =0.02 (fm" 3 ) N ^ 
1,1,1,1,1, 


- — 

' S 

" (b) \s, 

NX 

. p b =0.005 (fm" 3 ) 
1,1,1,1,1, 


*»* 

* **». t 

* N. — 

- (d?\N - 

. p b =0.02(fm 3 ) > 
1,1,1,1,1, 



10 



34567 34567 
T (MeV) 



FIG. 7: (Color online) The binding energy and free energy of 
40 Ca, produced in fragmentation of symmetric nuclear mat- 
ter. The left panels refer to pb = 0.005 fm -3 , the right panels 
correspond to pb = 0.02 fm~ 3 . The full (blue), dashed (red) 
and the dash-dot (black) lines have the same meaning as in 
Fig. 4. 



asymmetric nuclear matter at the two baryon densities 
we have considered for the three fragment species. We 
display only the full calculations (1,1,1). The blue full 
lines correspond to X = 0.0, the blue dotted lines re- 
fer to X = 0.2. 40 Ca and 56 Fe being symmetric and 
nearly symmetric nuclei, respectively have comparatively 
larger population in symmetric nuclear matter. For the 
more asymmetric 150 Sm nucleus, population is larger in 
the asymmetric matter. The difference between the two 
calculations at the two asymmetries is more pronounced 
at lower temperatures gradually narrowing down as the 
temperature is raised. The internal nucleon density of 
the fragments and their binding energies, however, show 
no significant change when the isospin asymmetry of the 
matter changes. This is displayed in Fig. 10 and Fig. 11, 
respectively for the three nuclear fragments at the two 
baryon densities; the blue full and dotted lines nearly 
overlap each other over the whole temperature range we 
work in. 
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FIG. 8: (Color online) Same as in Fig. 7 for 150 Sm. 
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FIG. 9: (Color online) Comparison of fragment population in 
symmetric and asymmetric nuclear matter for the case (1,1,1) 
calculated for two baryon densities. The full and dotted lines 
correspond to X — 0.0 and X = 0.2, respectively. 



IV. CONCLUDING REMARKS 

Analysis of nuclear multifragmentation data showed 
that for the results calculated in thermodynamic mod- 
els to conform to the experimentally observed ones, the 
established nuclear parameters taken as inputs in these 

1 _ j_* 1 _ _l _l lo-il lr*rj Inrl lr*/?fl n l_ 



a fact points out that the properties of the fragments 
produced in nuclear disassembly might have got modi- 
fied because of the interaction of the fragments with the 
embedding environment they are created in. The calcu- 
lations presented in this paper throw light in a quanti- 
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FIG. 10: (Color online) Comparison of the internal nucleon 
densities of the fragments produced in symmetric and asym- 
metric (X — 0.2) nuclear matter calculated for two baryon 
densities. The full and dotted lines have the same meaning 
as in Fig. 9. 
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gies of the fragments produced in symmetric and asymmetric 
(X = 0.2) nuclear matter calculated for two baryon densities. 
The full and dotted lines have the same meaning as in Fig. 9. 



be. For simplicity, the fragmented system was assumed 
to contain a collection of only one kind of nuclear species 
in addition to neutrons and protons. For comparison, 
a benchmark calculation of the hot fragments in phase- 
equilibrium with its own vapor was also done. In both 
cases, it was seen that the fragments expand with tem- 
perature as expected, but compared to the results from 
the phase-equilibrium (benchmark) calculation, the em- 
bedding medium in multifragmentation produced larger 
changes in the fragment properties. The fragments get 
comparably more stretched; the volume energy, surface 
properties or the symmetry properties of the fragments 
undergo the consequential changes. 

Questions may arise on the justification of the choice 
of only one kind of species in the calculations. Close 
examination of the results shows that in the medium, the 
interaction of the nuclcons with the fragments plays the 
dominant role in bringing forth the modification in the 
fragment properties. The fragment-fragment interaction 
has a very nominal role. The selection of the multispecies 
in the medium may not thus possibly alter the results 
much. 

Experiments with intermediate energy heavy ion 
beams in the last few decades have indicated [53[ that 
nuclei can sustain only temperatures that are much lower 
than the critical temperature (~ 16 MeV) for symmetric 
nuclear matter. The origin of such a limiting temperature 
is usually traced to an interplay between the Coulomb in- 
stability and the corrections due to the finite size of the 

nuclear drops. Phase-equilibrium calculations for hot iso- 
i i i i • l ii ii ■ _ •nil* ■ j ■ 



temperatures ~5— 7 MeV for heavier nuclei [Fll. IH2I |54| ■ 
In the calculations presented in this paper, it is seen that 
for nuclei dripped in a nuclear soup, the interaction with 
the surrounding medium might overcome the said insta- 
bility to a certain extent and extend somewhat the limit 
in temperature the nuclei may hold. This might be of sig- 
nificant relevance in the context of nuclear astrophysics 
and needs further exploration. 
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Appendix A 

From Eqs. (11) and (18), the density of the free nucle- 
ons p r is found to be 

Air , 
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Similarly, from Eqs. (12) and (19), the fragment density 
Pa comes out as 



p A = -^(2m A T) 3 / 2 J 1/2 ( VA ), 



or 



2tt 



p A = Z_ {2m * A T) 3 ^B 1/2 ( VA ), 



(A.2) 



(A.3) 



depending on whether the fragments are fermions or 
bosons. The Jk(v) an d Bk(rj) are the Fermi and Bose 
integrals, 



Bk(v) 



k dx 



e( x ~ r i) + 1 ' 
x k dx 



, e ( ;E -'')-l' 
The expression for V® in Eq. (22) is given as, 

V? = -47TO 3 {1 - d 2 {2 PN ) K ) (C lPr + C u p- T ) 

16w 2 a 3 



(A.4) 



b 2 h 3 

XJ3/2(V-t) 

An 



C i (2m*T) 5 / 2 J 3/2 (r ?T ) + C u {2m%Tf/ 2 



+ l^M^^ - l)(C lPlT + C uPl . T ) 



h 3 b 2 



[G(2m^T) 5 /V 3/2 (^) + C u (2mt\) 5/2 



XJsMvf-r)} + (ClPlr + C u pi- T )d 2 [p N + pi] K j. (A.5) 

In this equation, the fugacity 77^ is defined corresponding 
to the nucleons of density pi T inside the fragment (exactly 
in parallel definition of the fugacity r\ r of the free nucleons 
corresponding to the free nucleon density p T ); raf* is the 
effective mass of these nucleons. 

Similarly, expressions for V^, V 2 , V A and V A are, 

Vr = + C uP - r ] + \l{d + C u )^-, (A.6) 



1 \, v f An (2m*Tf/ 2 

T 

X-J3/2{Vr)+d 2 {p N +pi} K } 

+ j^{Ci(2mf T *Tf 2 J 3/2 ( V f T ) 
+C u (2mf* T Tf/ 2 J 3/2 (r ] f_ T )}' 
+IaPa'^2pit< Cipi T (-l 



< Pa> , 
b 2 b 2 b? 



PlT * V 



<Pa > 
b 2 



An {2mtTf/ 2 A 

W Wr 3/2{VlT > 



An {2mf* T Tf/< 



J3/2 {ritr) 



b 2 h? p,_ T 

+I A p A d 2 \ [ Y J Pir(Cipir + C u pi- T ) (2 Pl ) K , (A.8; 

r 

Vl = J2( C l^ + C uP i- T )(^- Plr + -jL T ). (A.9) 



In Eqs. (A.5)-(A.8), < p A > is the mean squared value 
of the fragment momentum in AN matter. Its value is 
given by, 



< p\ >= (2m* A T)C A , 



(A.10) 



where Ca is given by Eqs. (40) or (41). The integrals I 
and I a appearing in Eq. (A.5)-(A.9) for nucleon-nucleus 
and nucleus-nucleus interactions are given by, 



1=1 dr I dK 

Va 



e -|r+R|/a 

|r + R|/a' 



(A.ll) 



V 2 = Ana Z Kd 2 {2p N ) K - 1 Y,[CiPr'+C uP ^]p Tl 



+ -^IpaPn{CiPit + C u pi_ T ) 
x[Kd 2 ( PN + pi)"- 1 } , 



(A.7) 



p p „-|R+r-r'|/ a 

I A = I dr dr' / rfR— — , 

Iva Jva J |R + r-r'|/a 



(A.12) 



Integrations on R exclude the fragment volumes. The 
integrals are evaluated numerically. 
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